Effect of proline content and histidine ligation on the dynamics of Ω-loop D and the peroxidase activity of iso-1-cytochrome c.

To study how proline residues affect the dynamics of Ω-loop D (residues 70 to 85) of cytochrome c, we prepared G83P and G83A variants of yeast iso-1-cytochrome c (iso-1-Cytc) in the presence and absence of a K73H mutation. Ω-loop D is important in controlling both the electron transfer function of Cytc and the peroxidase activity of Cytc used in apoptosis because it provides the Met80 heme ligand. The G83P and G83A mutations have no effect on the global stability of iso-1-Cytc in presence or absence of the K73H mutation. However, both mutations destabilize the His73-mediated alkaline conformer relative to the native state. pH jump stopped-flow experiments show that the dynamics of the His73-mediated alkaline transition are significantly enhanced by the G83P mutation. Gated electron transfer studies show that the enhanced dynamics result from an increased rate of return to the native state, whereas the rate of loss of Met80 ligation is unchanged by the G83P mutation. Thus, the G83P substitution does not stiffen the conformation of the native state. Because bis-His heme ligation occurs when Cytc binds to cardiolipin-containing membranes, we studied the effect of His73 ligation on the peroxidase activity of Cytc, which acts as an early signal in apoptosis by causing oxygenation of cardiolipin. We find that the His73 alkaline conformer suppresses the peroxidase activity of Cytc. Thus, the bis-His ligated state of Cytc formed upon binding to cardiolipin is a negative effector for the peroxidase activity of Cytc early in apoptosis.

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST.

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

c-di-GMP regulates activity of the PlzA RNA chaperone from the Lyme disease spirochete.

PlzA is a c-di-GMP-binding protein crucial for adaptation of the Lyme disease spirochete Borrelia (Borreliella) burgdorferi during its enzootic life cycle. Unliganded apo-PlzA is important for vertebrate infection, while liganded holo-PlzA is important for survival in the tick; however, the biological function of PlzA has remained enigmatic. Here, we report that PlzA has RNA chaperone activity that is inhibited by c-di-GMP binding. Holo- and apo-PlzA bind RNA and accelerate RNA annealing, while only apo-PlzA can strand displace and unwind double-stranded RNA. Guided by the crystal structure of PlzA, we identified several key aromatic amino acids protruding from the N- and C-terminal domains that are required for RNA-binding and unwinding activity. Our findings illuminate c-di-GMP as a switch controlling the RNA chaperone activity of PlzA, and we propose that complex RNA-mediated modulatory mechanisms allow PlzA to regulate gene expression during both the vector and host phases of the B. burgdorferi life cycle.

Alkaline State of the Domain-Swapped Dimer of Human Cytochrome c: A Conformational Switch for Apoptotic Peroxidase Activity.

A 2.08 Å structure of an alkaline conformer of the domain-swapped dimer of K72A human cytochrome c (Cytc) crystallized at pH 9.9 is presented. In the structure, Lys79 is ligated to the heme. All other domain-swapped dimer structures of Cytc have water bound to this coordination site. Part of Ω-loop D (residues 70-85) forms a flexible linker between the subunits in other Cytc domain-swapped dimer structures but instead converts to a helix in the alkaline conformer of the dimer combining with the C-terminal helix to form two 26-residue helices that bracket both sides of the dimer. The alkaline transition of the K72A human dimer monitored at both 625 nm (high spin heme) and 695 nm (Met80 ligation) yields midpoint pH values of 6.6 and 7.6, respectively, showing that the Met80 → Lys79 and high spin to low spin transitions are distinct. The dimer peroxidase activity increases rapidly below pH 7, suggesting that population of the high spin form of the heme is what promotes peroxidase activity. Comparison of the structures of the alkaline dimer and the neutral pH dimer shows that the neutral pH conformer has a better electrostatic surface for binding to a cardiolipin-containing membrane and provides better access for small molecules to the heme iron. Given that the pH of mitochondrial cristae ranges from 6.9 to 7.2, the alkaline transition of the Cytc dimer could provide a conformational switch to tune the peroxidase activity of Cytc that oxygenates cardiolipin in the early stages of apoptosis.

Effect on intrinsic peroxidase activity of substituting coevolved residues from Ω-loop C of human cytochrome c into yeast iso-1-cytochrome c.

Naturally-occurring variants of human cytochrome c (Cytc) that induce thrombocytopenia IV occur within Ω-loop C (residues 40-57). These variants enhance the peroxidase activity of human Cytc apparently by facilitating access to the heme by destabilizing Ω-loops C and D (residues 70-85). Given the importance of peroxidase activity in the early stages of apoptosis, we identified three sites with the EVmutation algorithm in or near Ω-loop C that coevolve and differ between yeast iso-1-Cytc and human Cytc. We prepared iso-1-Cytc variants with all possible combinations of the S40T, V57I and N63T substitutions to determine if these residues decrease the peroxidase activity of iso-1-Cytc to that of human Cytc producing an effective off state for a peroxidase signaling switch. At pH 6 and above, all variants significantly decreased peroxidase activity. However, the correlation of peroxidase activity with local and global stability, expected if cooperative unfolding of Ω-loops C and D is required for peroxidase activity, was generally poor. The m-values derived from the guanidine hydrochloride dependence of the kinetics of imidazole binding to horse Cytc, which is well-characterized by native-state hydrogen exchange methods, and K72A/K73A/K79A iso-1-Cytc show that local structural fluctuations and not subglobal cooperative unfolding of Ω-loops C and D are sufficient to permit binding of a small molecule like peroxide to the heme. A 2.46 Å structure of N63T iso-1-Cytc identifies a change to a hydrogen bond network linking Ω-loops C and D that could modulate the local fluctuations needed for the intrinsic peroxidase activity of Cytc.

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A.

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔGu°'(H2O), of UBA(1) is 2.4 kcal mol-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, kf, in the absence of a denaturant of 13,000 s-1 and a Tanford ß-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone 15NH, 13CO, and 13Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues.

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

Curvature-Dependent Binding of Cytochrome c to Cardiolipin.

Cytochrome c binds cardiolipin on the concave surface of the inner mitochondrial membrane, before oxidizing the lipid and initiating the apoptotic pathway. This interaction has been studied in vitro, where mimicking the membrane curvature of the binding environment is difficult. Here we report binding to concave, cardiolipin-containing, membrane surfaces and compare findings to convex binding under the same conditions. For binding to the convex outer surface of cardiolipin-containing vesicles, a two-step structural rearrangement is observed with a small rearrangement detectable by Soret circular dichroism (CD) occurring at an exposed lipid-to-protein ratio (LPR) near 10 and partial unfolding detectable by Trp59 fluorescence occurring at an exposed LPR near 23. On the concave inner surface of cardiolipin-containing vesicles, the structural transitions monitored by Soret CD and Trp59 fluorescence are coincident and occur at an exposed LPR near 58. On the concave inner surface of mitochondrial cristae, we estimate the LPR of cardiolipin to cytochrome c is between 50 and 100. Thus, cytochrome c may have adapted to its native environment so that it can undergo a conformational change that switches on its peroxidase activity when it binds to CL-containing membranes in the cristae early in apoptosis. Our results show that membrane curvature qualitatively affects peripheral protein-lipid interactions and also highlights the disparity between in vitro binding studies and their physiological counterparts where cone-shaped lipids, like cardiolipin, are involved.

The Human Cytochrome c Domain-Swapped Dimer Facilitates Tight Regulation of Intrinsic Apoptosis.

Oxidation of cardiolipin (CL) by cytochrome c (cytc) has been proposed to initiate the intrinsic pathway of apoptosis. Domain-swapped dimer (DSD) conformations of cytc have been reported both by our laboratory and by others. The DSD is an alternate conformer of cytc that could oxygenate CL early in apoptosis. We demonstrate here that the cytc DSD has a set of properties that would provide tighter regulation of the intrinsic pathway. We show that the human DSD is kinetically more stable than horse and yeast DSDs. Circular dichroism data indicate that the DSD has a less asymmetric heme environment, similar to that seen when the monomeric protein binds to CL vesicles at high lipid:protein ratios. The dimer undergoes the alkaline conformational transition near pH 7.0, 2.5 pH units lower than that of the monomer. Data from fluorescence correlation spectroscopy and fluorescence anisotropy suggest that the alkaline transition of the DSD may act as a switch from a high affinity for CL nanodiscs at pH 7.4 to a much lower affinity at pH 8.0. Additionally, the peroxidase activity of the human DSD increases 7-fold compared to that of the monomer at pH 7 and 8, but by 14-fold at pH 6 when mixed Met80/H2O ligation replaces the lysine ligation of the alkaline state. We also present data that indicate that cytc binding shows a cooperative effect as the concentration of cytc is increased. The DSD appears to have evolved into a pH-inducible switch that provides a means to control activation of apoptosis near pH 7.0.

Cardiolipin Preferentially Partitions to the Inner Leaflet of Mixed Lipid Large Unilamellar Vesicles.

Cardiolipin (CL), an anionic phospholipid constituting 20% of the inner mitochondrial membrane (IMM) of eukaryotes, stabilizes electron transport chain (ETC) complexes and is a signaling agent in the early stages of apoptosis. For apoptosis, CL moves from the inner to the outer leaflet of the IMM via a poorly understood mechanism. Relative to cylindrically shaped lipids like dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylglycerol (DOPG), cone-shaped CL should prefer the concave surfaces of lipid bilayers. Using the fluorophore, 1,1,2,2-tetrakis[4-(2-trimethylammonioethoxy)phenyl]ethene, we have measured CL versus DOPG partitioning to the inner versus the outer leaflet of liposomes in mixed lipid systems with DOPC. DOPG shows no leaflet preference. However, CL has a 4:1 preference for the concave surface of the inner leaflet of liposomes. To further test the inner leaflet preference of CL, we show that cytochrome c binding to the outer convex surface of 20% CL/80% DOPC vesicles is strongly attenuated. Because the outer leaflet of intracristal regions of the IMM has a concave curvature, the preference of CL for concave surfaces may facilitate the transport of CL from the inner to the outer leaflet of the IMM needed for CL remodeling, the optimal functioning of the ETC, and signaling in the early stages of apoptosis.

Naturally Occurring A51V Variant of Human Cytochrome c Destabilizes the Native State and Enhances Peroxidase Activity.

The A51V variant of human cytochrome c is linked to thrombocytopenia 4 (THC4), a condition that causes decreased blood platelet counts. A 1.82 Å structure of the A51V variant shows only minor changes in tertiary structure relative to the wild-type (WT) protein. Guanidine hydrochloride denaturation demonstrates that the global stability of the A51V variant is 1.3 kcal/mol less than that of the WT protein. The midpoint pH, pH1/2, of the alkaline transition of the A51V variant is 1 unit less than that of the WT protein. Stopped-flow pH jump experiments show that the A51V substitution affects the triggering ionization for one of two kinetically distinguishable alkaline conformers and enhances the accessibility of a high-spin heme transient. The pH1/2 for acid unfolding of the A51V variant is 0.7 units higher than for that of the WT protein. Consistent with the greater accessibility of non-native conformers for the A51V variant, the kcat values for its peroxidase activity increase by 6- to 15-fold in the pH range of 5-8 versus those of the WT protein. These data along with previously reported data for the other THC4-linked variants, G41S and Y48H, underscore the role of Ω-loop C (residues 40-57) in modulating the peroxidase activity of cytochrome c early in apoptosis.

Effect of V83G and I81A Substitutions to Human Cytochrome c on Acid Unfolding and Peroxidase Activity below a Neutral pH.

Mitochondrial cytochrome c is a highly conserved protein in eukaryotes. Certain functions of cytochrome c have been tuned during evolution. For instance, the intrinsic peroxidase activity of human cytochrome c is much lower than that of the yeast counterpart. Structural studies on K72A yeast iso-1-cytochrome c [McClelland, L. J., et al. (2014) Proc. Natl. Acad. Sci. USA, 111, 6648-6653] revealed that residues 81 and 83 in Ω-loop D (residues 70-85) may be gatekeeper residues for the peroxidase activity linked to intrinsic apoptosis. Amino acids at both positions evolve to more sterically demanding amino acids. We hypothesized that residues 81 and 83 evolved such that steric constraints at these positions tune down the peroxidase activity of human cytochrome c. To test this hypothesis, I81A and V83G variants of human cytochrome c were prepared. Our results show that the I81A substitution significantly influences the thermodynamics and kinetics of access to alternate conformers of human cytochrome c, while the V83G substitution has a more modest effect on these properties. The I81A variant also shows a significant enhancement in peroxidase activity, particularly below pH 7, whereas the V83G variant has a similar peroxidase activity to wild-type human cytochrome c. However, neither variant increases the peroxidase activity of human cytochrome c to the level of yeast iso-1-cytochrome c, indicating that other substructures of cytochrome c are also involved in tuning the intrinsic peroxidase activity of mitochondrial cytochrome c.

Electrostatic Constituents of the Interaction of Cardiolipin with Site A of Cytochrome c.

Cytochrome c binds to cardiolipin (CL) on the inner mitochondrial membrane during the initial stages of apoptosis where it oxidizes CL, promoting its release into the cytoplasm where it initiates apoptosis. Previous work has identified interaction sites on cytochrome c involved in the cytochrome c-CL interaction. The contributions of the lysines attributed to site A, the anionic site, are studied here to elucidate the relative importance of each for electrostatic interaction of cytochrome c with CL at pH 8, conditions where site A is dominant. A set of single, double, and quadruple lysine to alanine variants of yeast iso-1-cytochrome c, at sequence positions 72, 73, 86, and 87, show that all contribute to the site A-mediated interaction with CL. All variants experience two sequential structural rearrangements as the lipid to protein ratio (LPR) increases. At a low LPR near 10, all variants undergo a small heme-centered structural change detected by Soret circular dichroism. At higher LPRs ranging from 22 to 34, all variants partially unfold as detected by Trp59 emission. The robustness of the mechanism of interaction to sequential neutralization of the four lysines assigned to site A demonstrates that site A is more extensive than previously supposed. The nature of both structural rearrangements also depends on which lysines constitute site A. The peroxidase activity of cytochrome c in the early stages of apoptosis depends on the nature of structural rearrangement near the heme. Thus, the lysines that comprise site A may have evolved to optimize the peroxidase signaling switch.

Rapid quantification of vesicle concentration for DOPG/DOPC and Cardiolipin/DOPC mixed lipid systems of variable composition.

A novel approach to quantify mixed lipid systems is described. Traditional approaches to lipid vesicle quantification are time consuming, require large amounts of material and are destructive. We extend our recently described method for quantification of pure lipid systems to mixed lipid systems. The method only requires a UV-Vis spectrometer and does not destroy sample. Mie scattering data from absorbance measurements are used as input into a Matlab program to calculate the total vesicle concentration and the concentrations of each lipid in the mixed lipid system. The technique is fast and accurate, which is essential for analytical lipid binding experiments.

Humanlike substitutions to Ω-loop D of yeast iso-1-cytochrome c only modestly affect dynamics and peroxidase activity.

Structural studies of yeast iso-1-cytochrome c (L.J. McClelland, T.-C. Mou, M.E. Jeakins-Cooley, S.R. Sprang, B.E. Bowler, Proc. Natl. Acad. Sci. U.S.A. 111 (2014) 6648-6653) show that modest movement of Ω-loop D (residues 70-85, average RMSD versus the native structure: 0.81 Å) permits loss of Met80-heme ligation creating an available coordination site to catalyze the peroxidase activity mediated by cytochrome c early in apoptosis. However, Ala81 and Gly83 move significantly (RMSDs of 2.18 and 1.26 Å, respectively). Ala81 and Gly83 evolve to Ile and Val, respectively, in human cytochrome c and peroxidase activity decreases 25-fold relative to the yeast protein at pH 7. To test the hypothesis that these residues evolved to restrict the peroxidase activity of cytochrome c, A81I and G83V variants of yeast iso-1-cytochrome c were prepared. For both variants, the apparent pKa of the alkaline transition increases by 0.2 to 0.3 relative to the wild type (WT) protein and the rate of opening the heme crevice is slowed. The cooperativity of acid unfolding is decreased for the G83V variant. At pH 7 and 8, the catalytic rate constant, kcat, for the peroxidase activity of both variants decreases relative to WT, consistent with the effects on alkaline isomerization. Below pH 7, the loss in the cooperativity of acid unfolding causes kcat for peroxidase activity to increase for the G83V variant relative to WT. Neither variant decreases kcat to the level of the human protein, indicating that other residues also contribute to the low peroxidase activity of human cytochrome c.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'.

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein.

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c.

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign 1H and 15N backbone amide and 13C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Site A-Mediated Partial Unfolding of Cytochrome c on Cardiolipin Vesicles Is Species-Dependent and Does Not Require Lys72.

Measurements at pH 8 allow evaluation of binding of 100% cardiolipin vesicles to site A of cytochrome c without interference from other known binding sites. Site A encompasses Lys72, Lys73, Lys86, and Lys87, located in or adjacent to Ω-loop D (residues 70-85), which positions Met80 for binding to the heme. Binding of cytochrome c to cardiolipin disrupts Met80 heme binding, permitting peroxidase activity. Binding of cardiolipin to yeast iso-1-cytochrome c versus human cytochrome c is compared to assess how binding of cardiolipin to site A has evolved for cytochrome c from species that do not have a complete intrinsic apoptotic pathway to species that do. Using a nondestructive method of quantifying cardiolipin concentration, highly reproducible binding curves are obtained. The results indicate two sequential structural rearrangements on the surface of 100% cardiolipin vesicles. The first, more modest, structural rearrangement occurs at an exposed (outer leaflet) lipid:protein ratio of 8-10 for both cytochromes c. The second, occurring at higher lipid:protein ratios, causes significant unfolding of cytochrome c and requires a much higher lipid:protein ratio for human versus yeast cytochrome c. Higher lipid:protein ratios enhance the peroxidase activity of cytochrome c, suggesting that human cytochrome c has evolved a more stringent on/off switch for cardiolipin peroxidation in the early stages of apoptosis. For both human and yeast cytochrome c, the K72A mutation has only minor effects on binding to site A, suggesting that other nearby lysines can compensate for the lack of Lys72.

Effect of a K72A Mutation on the Structure, Stability, Dynamics, and Peroxidase Activity of Human Cytochrome c.

We test the hypothesis that Lys72 suppresses the intrinsic peroxidase activity of human cytochrome c, as observed previously for yeast iso-1-cytochrome c [McClelland, L. J., et al. (2014) Proc. Natl. Acad. Sci. U. S. A. 111, 6648-6653]. A 1.25 Å X-ray structure of K72A human cytochrome c shows that the mutation minimally affects structure. Guanidine hydrochloride denaturation demonstrates that the K72A mutation increases global stability by 0.5 kcal/mol. The K72A mutation also increases the apparent pKa of the alkaline transition, a measure of the stability of the heme crevice, by 0.5 unit. Consistent with the increase in the apparent pKa, the rate of formation of the dominant alkaline conformer decreases, and this conformer is no longer stabilized by proline isomerization. Peroxidase activity measurements show that the K72A mutation increases kcat by 1.6-4-fold at pH 7-10, an effect larger than that seen for the yeast protein. X-ray structures of wild type and K72A human cytochrome c indicate that direct interactions of Lys72 with the far side of Ω-loop D, which are seen in X-ray structures of horse and yeast cytochrome c and could suppress peroxidase activity, are lacking. Instead, we propose that the stronger effect of the K72A mutation on the peroxidase activity of human versus yeast cytochrome c results from relief of steric interactions between the side chains at positions 72 and 81 (Ile in human vs Ala in yeast), which suppress the dynamics of Ω-loop D necessary for the intrinsic peroxidase activity of cytochrome c.